Large- J approach to strongly coupled spin-orbital systems 
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We present a novel approach to study the ground state and elementary excitations in compounds 
where spins and orbitals are entangled by on-site relativistic spin-orbit interaction. The appropriate 
degrees of freedom are localized states with an effective angular momentum J. We generalize J 
to arbitrary large values while maintaining the delicate spin-orbital entanglement. After projecting 
the inter-site exchange interaction to the manifold of effective spins, a systematic 1/ J expansion of 
the effective Hamiltonian is realized using the Holstein-Primakoff transformation. Applications to 
representative compounds Sr2ir04 and particularly vanadium spinels AV2O4 are discussed. 



Transition metal compounds with partially filled or- 
bitals have attracted considerable attention in recent 
years^ The interplay of spin and orbital degrees of free- 
dom leads to a variety of interesting ground states and 
elementary excitations. Particularly, a long-range orbital 
order in the ground state makes possible the propagating 
orbital excitations, or orbitons^ in addition to the famil- 
iar spin waves. The intricate spin-orbital interaction dic- 
tated by the Kugel-Khomskii 3 - type superexchange (SE) 
Hamiltonian gives rise to a magnon spectrum which de- 
pends on the orbital order, and vice versa. Meanwhile, 
interaction between the two types of quasiparticles leads 
to their scattering and decay. 

Recently, there is a growing number of orbitally degen- 
erate compounds whose elementary excitations can not 
be described by the above paradigm based on magnons 
and orbitons. A common feature shared by these com- 
pounds is the presence of a strong relativistic spin-orbit 
(SO) interaction Vls — A(L ' S); most notable among 
them are vanadium and iridium oxides containing V 3+ 
and Ir 4+ ions, respectively4iSiSil To understand these 
systems, one should start with localized degrees of free- 
dom which diagonalize the SO coupling^^S Such atomic 
states are usually composed of complex orbitals and are 
characterized by an effective angular momentum J. The 
on-site spin-orbital entanglement leads not only to novel 
ground states with orbital-moment orderingj£ii£ but also 
to the appearance of new types of quasiparticles which 
carry both spin and orbital characters. 

Although several theoretical studies have been devoted 
to understanding the ground-state properties in the limit 
of strong SO coupling 1 1 1 1 1 1 2 very few works address 
the problem of elementary excitations. A straightfor- 
ward approach is to use the magnon/orbiton description 
starting from the A = ground state and to treat SO 
interaction as a perturbation^ However, as the large- 
A ground state is fundamentally different from that of 
A = 0, such a formulation would be difficult to capture 
features like ordering of orbital moments. Another alter- 
native is the magnetic-exciton model. Its application to 
vanadium spinels reveals an excitation spectrum charac- 
terized by transitions between states with different J z .— 
Despite its success in providing a basic structure of exci- 
tations in the large A limit, the magnetic-exciton model 



relies on an ill-controlled mean-field decoupling of the SE 
interaction. Besides, the model itself does not offer much 
insight to the ground-state structure. 

In this paper, we propose a new theoretical framework 
to study magnetically ordered ground state and its quasi- 
particles in systems with a dominating SO interaction. 
Instead of naively extending the length of the effective 
spin, our large- J generalization preserves the delicate en- 
tanglement between spins and orbitals. After obtaining 
the classical ground states in the J — * 00 limit, a sys- 
tematic 1/J expansion can be attained with the aid of 
Holstein-Primakoff transformation. Using the large- J ap- 
proach, we examine two canonical examples: the iridium 
and vanadium compounds with effective spins J = 1/2 
and J = 2, respectively. In the case of Ir 4+ ion, the large- 
J approach is equivalent to conventional large-5* method; 
the classical Hamiltonian has the same form as its quan- 
tum counterpart. 

The situation for V 3+ ion is more complicated. By 
minimizing the J — > 00 limit of the effective Hamilto- 
nian, we obtain a classical phase diagram of vanadium 
spinels. Our result provides a coherent explanation for 
the collinear and orthogonal antifcrromagnetic orders ob- 
served in spinels ZnV2 04 and MnV204, respectively. 
The SE interaction in vanadium spinels has rather dis- 
tinct characteristics in spin and orbital channels: while 
the spin exchange preserves a SU(2) symmetry, the or- 
bital part resembles an anisotropic 3-state Potts model. 
We find that the strong spatial dependence of the orbital 
exchange leads to a gapped excitation spectrum with a 
narrow bandwidth. 

Spin-orbit coupling and effective spins. We start by 
describing the general large-J extension and then dis- 
cuss specific applications to Ir 4+ and V 3+ ions with 
d 5 and d 2 configurations, respectively. In both cases, 
the orbital degrees of freedom of a partially filled ti a 
triplet are described by an effective angular momentum 
of length L' = 1. The true angular momentum is given 
by L = cvL', where a of order one is the so-called co- 
valency factor^ The SO interaction Vls = ce\(L' ■ S) 
is diagonalized by the eigenstates of the "total angular 
momentum" J = L + S. Depending on the sign of a, 
the ground state has an effective angular momentum of 
length J ff = L' ± S and is separated from the excited 
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states by a gap of order A. 

In order to carry out a controlled I/J expansion of 
the effective Hamiltonian, we first generalize J to arbi- 
trary values by considering a fictitious ion with N iden- 
tical copies of t 2g triplets. The generalized SO interac- 
tion is Vls = oz\J2 m =i ' S m , where L' m and S m are 
orbital momentum and spin operators of the m-th t 2g 
triplet, respectively. We further assume that only states 
which are symmetric with respect to the N degenerate 
t 2g triplets are allowed physical states. In the atomic 
ground state, each t 2g triplet has an angular momentum 
of length J m — J c ff = V ± S, where 3 m = \J m + S m , and 

the symmetric sum J = X)m=i nas ^ ne maximum 
length J = NJ c g. 

Single-ion operators such as orbital projections and 
spins are extended accordingly: P a = 53m=i Pm,a an d 
S = 5Zm=i To obtain their representation in the 

J = NJ c g subspace, we expand its basis | J, J z ) with re- 
spect to the quantum number J mjZ of the m-th tig triplet 

\J,Jz) = ^2Cj^ Jm J Jeff, J m , z )®| J" J c ff, J z - J m , z ),(l) 



where Cj Z! j m z are the Clebsch-Gordon coefficients. By 
further expressing | J m , J mjZ ) in terms of \L' m z , 

Sm,z} i one 

obtains the matrix elements of operators acting on the fri- 
th tig orbitals. As the J = NJ e g subspace only contains 
permutationally symmetric states, the overall operator 
at a site is just N times that of a single t 2g triplet, e.g. 
S = NS m . 

Taking into account the inter-site exchange interaction, 
the effective Hamiltonian of the angular momenta J^ in 
the A — ► oo limit is given by H c g = VHseP, where V is 
the projection operator of the J = NJ c g manifold, and 
-Hse is the spin-orbital superexchange Hamiltonian. 

7r 4+ ions: L' = 1, S = 1/2. The iridium ions in 
compounds such as Sr2lr04 are in the low spin 5d 5 con- 
figuration with a positive covalency a « If** the atomic 
ground state has an effective spin J c ff = 1/2. By in- 
troducing N identical t 2g replicas as described above, 
we generalize the effective spin to arbitrary large values 
J = N/2. Restricted to the lowest-energy manifold, we 
find P xy — P yz = P zx = 2J/3, i.e. the electrons equally 
populate the three t 2g orbitals. The angular momentum 
and spin operators are given by 



-J/3, 



L' = 4J/3. 



(2) 



The magnetic moment of the Ir 4+ ion is /i = ^s(2S + 
aL 1 ) w -^p-J. For compounds with 180° Ir-O-Ir bonds, 
the projected Hamiltonian 

ffeff = Jx J i ■ J J ( Ji ' ^ ) ( J J ' f « ) ( 3 ) 

is dominated by a Heisenberg isotropic exchange plus a 
pesudo-dipolar interaction^ Here J\ w AJ and J 2 w 
277,/ are effective exchange constants, 77 = Jh/U is the 



ratio of Hund's coupling to on-site Coulomb repulsion U, 
and J — t 2 /U defines the overall exchange energy scale. 
More importantly, we find that the large- J version of the 
effective Hamiltonian, and particularly the classical limit 
(J — * 00), has the same form as the quantum Hamilto- 
nian first derived in Ref. [9|. The large- J approach in 
this case thus is equivalent to the conventional large-S" 
extension. In Sr2ir04, the 'weak' ferromagnetic moment 
accompanying the ground-state antiferromagnetic order 
is explained by treating the spins as classical objects and 
taking into account the staggered rotations of the IrC>6 
octahedral The large- J approach thus provides a the- 
oretical basis for the classical treatment of the effective 
quantum Hamiltonian. 

V 3+ ions: L' = 1, S = 1. Vanadium ions in cubic 
and spinel vanadates have a 3d 2 configuration and a neg- 
ative a. Consequently, the effective angular momentum 
minimizing the SO interaction has a length of J c g — 2. 
The large- J generalization leads to the following orbital 
projection operator 



_ J(J-1) + J Z 2 
xy 2(2J-1) ' 



(4) 



For projections P yz / zx , we replace J z by J x and J y , re- 
spectively. It could be easily checked that P yz + P zx + 
Pxy = J, which is just the total number of electrons 27V. 
Similarly, projected to the J = 2N subspace, the spin 
and orbital momentum operators are 



J/2, L' = J/2. 



(5) 



The angular momentum quanta are evenly divided be- 
tween the spin and orbital channels. The V 3+ ion has a 
reduced magnetic moment fi — ^s(l — |a|/2)J « ^J- 
Due to the spin-orbital entanglement of the J = 2N 
states, the representation of the operator product P a S 
in this subspace is different from the matrix product of 
individual operators. Instead, we find 



PxyS 



Js J J z 



(J-l)J 

2(2J-1) 2(J-1)(2J-1)' 



(6) 



The cubic symmetry of the system allows us to obtain 
the expressions with P yz and P zx projections by simply 
replacing J z by J x and J y , respectively. 

The discussion so far has been general and is applica- 
ble to vanadium compounds with V 3+ ions in an octa- 
hedral crystal field. We now focus on the case of spinels 
where the vanadium ions form the pyrochlore lattice, a 
network of corner-sharing tetrahedra. The 90° V-O-V 
bonds make direct exchange the primary mechanism of 
the SE interactionj 11 ' 16 The resulting Potts-like orbital 
interaction leads to a highly anisotropic effective Hamil- 
tonian in both coordinate and spin spaces. 

Restricted to the J = 27V manifold, the effective quan- 
tum Hamiltonian can be obtained by substituting the 
operator expressions (HI)-© into the SE Hamiltonian for 
vanadium spinels. After introducing normalized classical 
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spins Ji = Jj/ J, the effective Hamiltonian in the classicai 
limit J — > oo becomes 

<«> 

+ l 7 2 J 2 ^(l + J^ 7 )(l + J| 7 )J i -J j (7) 

(«> 

— JzJ 2, [l — («A,7^J,7) 2 ] Jj ' Jj- 

Here the subscript 7 = 7(1,7) denotes the x, y ot z compo- 
nent of the angular momentum Jj depending on whether 
the nearest-neighbor bond (ij) points along the (Oil), 
(101), or (110) directions, respectively. To lowest order in 
T), the various exchange constants are J\ ~ ^(1 + 5r))J : 
J 2 ~ i(l - 77)J, and J 3 - I 

Classical ground state of vanadium spinels. There has 
been much experimental effort to understand the proper- 
ties of spinel vanadates ^4V2C>4. For divalent ion A = Zn, 
Cd, and Mg, magnetic moments of V 3+ ions are ordered 
collinearly in the ground stated Theoretical models 
based on a large SO coupling have been proposed to ex- 
plain the observed spin and orbital orderingi&iiiiS How- 
ever, a recent experimental characterization of another 
spinel MnV204 shows an orthogonal magnetic structure 
of vanadium spins £ The mechanisms that stabilize the 
orthogonal magnetic order remains unclear. Based on 
the large- J approach, here we propose a unifying model 
for the magnetic and orbital order in vanadium spinels. 

The partially filled ti g orbitals also couple to the 
distortions of the surrounding VOg octahedron. In- 
deed, both types of vanadium spinels undergo a cubic-to- 
tetragonal structural transition which accompanies the 
orbital ordering. The lattice distortion results in a split- 
ting of the tig triplet: Vjt = ^pSpPp, where j3 runs 
over the three tig states. Using Eq. Q, the JT coupling 
can be recast into 

Vjt = 5 X (J 2 + J 2 - 2J 2 )/V6 + 5 2 (J 2 X - J 2 )/V2, (8) 

where (5%, 62) transforming as a doublet irreducible rep- 
resentation under the symmetry group Oh denote the 
tetragonal and orthorhombic distortions of the VOq oc- 
tahedron, respectively. 

Structural distortions compatible with the observed 
tetragonal symmetry (lattice constants a = b > c) 
in general involves q = phonons with symmetries 
Eg and Fi g iH The E g distortion gives rise to a uni- 
form anisotropy S±, whereas softening of rotational F\ g 
phonons creates a staggered orthorhombic distortions 
±82 on VC>6 octahedra along the [110] and [110] chains, 
respectively. 17 Classical ground state of a tetrahedron, 
the building block of the pyrochlore lattice, is obtained 
by minimizing the sum of the effective Hamiltonian ([7]) 
and JT coupling |8]). The resulting phase diagram is 
shown in Fig. [Ha). 

The orthogonal magnetic order is stabilized at large 
values of the staggered F\ g distortion 82- Interestingly, 




FIG. 1: (a) Classical phase diagram of a tetrahedron. Black 
and white circles denote ±J Z components, respectively. Si 
and 82 denote the uniform E g and staggered F\ g distor- 
tions, respectively. A narrow low-symmetric phase with non- 
coplanar spins exists between phases B and C. (b) The q = 
orthogonal order and (c) the q = (001) collinear order are 
observed in spinels MnV^CU and ZnVaO.4, respectively. 



the orthogonal structure also minimizes the J\ term of 
if eg, which originates from the antiferro-orbital interac- 
tion in the original SE Hamiltonian. From Eq. ([3]), we 
see that d zx and d yz orbitals are fully occupied along the 
[110] and [110] chains, respectively. This antiferro-orbital 
ordering is consistent with the low-temperature symme- 
try IAi/a of MnV 2 4 ^ We note that the inclusion of SO 
interaction does not seem to affect the orbital-ordering 
pattern^ This could be attributed to the already mixed 
orbital states favored by a strong local trigonal distortion 
in MnV 2 4 . 

On the other hand, the collinear antifcrromagnetic or- 
der becomes the ground state when lattice distortion is 
dominated by the E g mode with <5i > 0, which gives rise 
to an easy-axis anisotropy according to ((SJ . The J2 term 
in ([7]) indicates that collinear spins parallel to the z-axis 
have the largest antiferromagnetic coupling on nearest- 
neighbor bonds lying in the xy planes. The ground state 
can then be viewed as a collection of antiferromagnetic 
Ising chains running along the [110] and [110] directions 
[Fig. [He)], as was indeed observed in ZnV^O^ The re- 
sulting orbital ordering has occupation numbers n xy — 1 
and n yz = n zx = 1/2, also consistent with the observed 
IAi/amd symmetry in ZhV^O^ A similar staggered or- 
dering of effective moment J is found in cubic vanadates 
at large X^- 

Holstein-Primakoff transformation. Once the classical 
ground state is determined from the effective Hamilto- 
nian, one can carry out a controlled 1/J expansion using 
the Holstein-Primakoff transformation: 

J z = J — a f a, J+ = \f2J~~aYaa, J_ = j\. (9) 
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FIG. 2: Quasiparticle spectra of vanadium spinels with (a) 
q = orthogonal and (b) q = (001) collinear magnetic orders 
in the limit of strong spin-orbit coupling (solid lines). Also 
shown for comparison is the magnon spectrum of the collinear 
order without spin-orbit coupling (dashed lines). In obtaining 
these spectra, we have used Si = 0, 82 = 0.2J and Si = 82 = 
in (a) and (b), respectively. The symmetry points in k- 
space are T = (0,0,0), L = (±, §, ±), Y = (0, 1,0), and Y' = 
(0, 1, |). The energy cuk is measured in units of J = t 2 /U. 
The inset in (a) shows the energy gap A as a function of 
staggered distortion 82- 

The linear quasiparticle spectrum can then be obtained 
with the aid of Bogoliubov transformation; the calculated 
spectra for the orthogonal and collinear magnetic orders 
are shown by the solid lines in Fig. [2j 

In the case of q = orthogonal magnetic order, we 
obtain an acoustic-like band whose gap scales as A a 
^82 — 5 2c in the vicinity of <5 2c = 0.154*7. The gap- 
less mode at 82c signals the transition from phases D 
to C in Fig. [1Tb). Although a similar gapped acoustic 
mode (A ss 1.5 meV) whose origin is attributed to the 
orthogonal order was indeed observed in MnV^O,^ the 
effect of the magnetic Mn 2+ ions (S = 5/2) remains to be 
clarified. On the contrary, all quasiparticle bands of the 
collinear magnetic order have an energy gap A ss 0.6J7 
even when Si = 0. This is because the staggered arrange- 
ment of frustrated bonds (parallel spins) in the q = (001) 
collinear order prevents the onset of the soft mode B 
shown in Fig. [T] 

Also shown for comparison is magnon spectrum of the 
collinear order in the absence of SO coupling (dashed 
lines) ii£ The Goldstone mode at the zone center reflects 



the global 0(3) symmetry of the superexchange Hamil- 
tonian in the A = limit. The rather narrow bandwidth 
of the excitation spectrum in the large A limit can be 
attributed to the fact that quasiparticles also carry or- 
bital degrees of freedom whose Potts-like interactions are 
static. 

Conclusion and outlook. To summarize, we present a 
new approach to the non-trivial problem of elementary 
excitations in systems with on-site spin-orbital entangle- 
ment. For such compounds, the ground-state magnetic 
order is usually accompanied by a long-range ordering 
of orbital moments. The bosonic elementary excitations 
resulting from deviations of the perfect order thus carry 
both spin and orbital degrees of freedom. Starting with 
the J — ► 00 limit, a systematic treatment of the elemen- 
tary excitations can be realized through the standard 1/ J 
expansion. The quasiparticle spectrum is obtained by di- 
agonalizing the quadratic Hamiltonian. Higher-order ef- 
fects such as quasiparticle scattering can also be studied 
using well developed methods from conventional large-5 
approach. 

It is worth noting that spin-orbital entanglement could 
exist even in the absence of SO coupling due to the opera- 
tor structure of SE Hamiltonian. For example, it is shown 
that composite spin-orbital fluctuations in cubic vana- 
dates lead to orbital-disordered phase and violation of the 
Goodenough-Kanamori rules. — As the large- J Hamilto- 
nian H e g = VHseP represents the lowest-order term in 
a J /A expansion, another important front is to include 
the effects of finite A. Higher-order terms can be ob- 
tained following the perturbation expansion of the SE 
Hamiltonian. In studying real compounds, A is usually 
of the same order of exchange J . Whether the excita- 
tions can be described by the large- J approach depends 
on the nature of the ground state. In the presence of 
orbital-moment ordering, we think the large-J method 
would be the appropriate choice. 

Finally, we note that the classical effective Hamilto- 
nian obtained in the J — > 00 limit provides a practical 
working model for analyzing the experimental data. Fur- 
thermore, being a generic approach, our method can be 
easily adapted to other compounds where a large SO cou- 
pling dominates the low-energy physics. 
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